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ABSTRACT 

We present detailed observations of a 2; ~ 1.99 cluster of submillimeter galaxies (SMGs), discovered 
as the strongest redshift spike in our entire survey of ^100 SMGs across 800 square arcmin. It is 
the largest blank-field SMG concentration currently known and has < 0.01% chance of being drawn 
from the underlying selection function for SMGs. We have compared UV observations of galaxies at 
this redshift, where we find a much less dramatic overdensity, having an 11% chance of being drawn 
from its selection function. We use this z ~ 1.99 overdensity to compare the biasing of UV- and 
submm-sclcctcd galaxies, and test whether SMGs could reside in less ovcrdcnsc environments, with 
their apparent clustering signal being dominated by highly active merger periods in modest mass 
structures. We discuss the probable mechanisms for the apparently different bias we see at the two 
wavelengths. This impressively active formation phase in a low mass cluster is not something seen 
in simulations, although we propose a toy model using merger bias which could account for the bias 
seen in the SMGs. While enhanced buildup of stellar mass appears characteristic of other high— z 
galaxy clusters, neither the UV- nor submm-galaxies in this structure exhibit larger stellar masses 
than their field galaxy counterparts (although the excess of SMGs in the structure represents a larger 
volume-averaged stellar mass than the field). Our findings have strong implications for future surveys 
for high-z galaxies at long wavelengths such as SCUBA2 and Herschel. We suggest that since these 
surveys will select galaxies during their episodes of peak starbursts, they could probe a much wider 
range of environments than just the progenitors of rich clusters, revealing more completely the key 
events and stages in galaxy formation and assembly. 
Subject headings: galaxies: optically faint radio — galaxies: evolution 



1. INTRODUCTION 

Submillimeter galaxies (SMGs) have provided an ef- 
ficient probe of ultraluminous (> 10"'^^ Lq) star forma- 
tion activity in the distant Universe (e.g, Smail, Ivison & 
Blain 1997: Blain et al. 2002; Borys et al. 2003; Coppin 
et al. 2008; Knudsen et al. 2008), with the bright submm 
detection (rest frame far-infrared - FIR) providing an 
unambiguous signal of large dust masses heated predom- 
inantly by young stars rather than AGN (e.g., Almaini 
et al. 1999; Fabian et al. 2000; Alexander et al. 2005; 
Menendez-Delmetre et al. 2007; Pope et al. 2008). Piec- 
ing together the detailed properties of SMGs has been 
difiicult because of the inherent dust obscuration to the 
regions of the largest bolometric output. High spatial 
resolution radio observations using the MERLIN inter- 
ferometer (Chapman et al. 2004a; Biggs & Ivison 2008) 
provided early indications that the dust and gas in SMGs 
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was more compact than typical galaxies, yet still resolved 
on scales three times larger than in local ULIRGs. Once 
verified directly through high resolution CO gas studies 
of SMGs (Tacconi et al. 2006, 2008) this was used to ar- 
gue that SMGs are scaled up ULIRGs. This radio and 
CO mapping of SMGs has suggested that the locus of 
far-IR emission is often offset from UV/optical emission, 
consistent with findings (Chapman et al. 2005 - COS) 
that even dust-corrected UV estimates of star formation 
rates in SMGs under-predict the true star formation rates 
(SFRs) by factors > 100 x on average. CO linewidths 
for a sizable sample of SMGs (Fraycr et al. 1998, Neri 
et al. 2003, Grcve et al. 2005) have suggested large dy- 
namical masses, > IO^^Mq, though not nearly as large 
as the halos in which they reside based on their cluster- 
ing properties (Blain et al. 2004). The large dynamical 
masses are associated with large (1 — 2 x 10"'^"'^ M0) stel- 
lar masses (Borys et al. 2005; Hainline 2008; Swinbank 
et al. 2008). Finally, the gas masses and SFRs define an 
average gas exhaustion timescale <50 Myrs, suggesting a 
duty cycle of ~10 for an underlying population of galax- 
ies hosting the submm-luminous events. These studies 
bring together a picture of SMGs as typically represent- 
ing massive, gas rich galaxies, which ought to be building 
elliptical galaxies that are found preferentially in clusters 
at the present day (Swinbank et al. 2006, 2008). 

The strongly clustered populations of red X-selccted 
galaxies at 2 = 2 — 4 (e.g., Daddi et al. 2003) and of lu- 
minous high- 2; radio galaxies (HzRG - Brand et al. 2005) 
suggest they are formed in the highest density perturba- 
tions at early epochs, implying that they are the progen- 
itors of local massive early-type galaxies and z ^ 1 ex- 
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Fig. 1. — Positions of all UV-selected galaxies an SMGs/SFRGs 
in the survey region are shown as small symbols, while those lying 
in the z = 1.99 structure are large symbols. Offset projected Mpc 
are shown on the axes for reference at z = 2. Galaxies in the z = 
1.99 structure is generally widely separated on the sky, with >Mpc 
separations between most SMGs/SFRGs except for two close (~ 
100 kpc) SMG/SFRG pairs in two cases. 



tremely red objects (EROs) which themselves have sim- 
ilarly strong clustering (e.g., Brown et al. 2005). HzRGs 
are believed to host massive black holes; thus may rep- 
resent some of the most massive galaxies at these epochs 
(e.g., de Breuck et al. 2000) and could trace some of the 
highest density environments. The identication of mod- 
est excesses of companion galaxies around HzRGs in the 
optical/near-infrared (e.g., Rottgering et al. 1996), X-ray 
(Pentericci et al. 2002; Small et al. 2003a), and submm 
(Ivison et al. 2000; Stevens et al. 2003), supports their 
association with overdense structures. 

Studies to date have suggested that the most FIR- 
luminous, distant galaxies are also clustered very 
strongly (Blain et al. 2004; Borys et al. 2004; Scott et al. 
2002, 2006; van Kampen et al. 2005; Farrah et al. 2006; 
Blake et al. 2006; GiUi et al. 2007). At face value this 
would indicate that they occupy the largest dark-matter 
halos ('^10^^ M0). However selection of galaxies by the 
most FIR-luminous specimens may be prone to finding 
the brief periods when systems of lower matter overden- 
sity are rapidly evolving (e.g. Scannapieco & Thacker 
2003; Furlanetto & Kamionkowski 2006). Indeed burst 
timcscales as short as 10 Myr have been estimated for 
some SMGs (e.g., Small et al. 2003b). Blain et al. (2004) 
calculated that, complex biases aside, SMGs have a clus- 
tering length which is consistent with a form of evolution 
ensuring their properties subsequently match the cluster- 
ing length typical of evolved red galaxies at 2 1 and 
finally of rich clusters of galaxies at the present epoch. 

However, as pointed out by Blain et al. (2004), this is 
inconsistent with the expected density of rich clusters. 
It is very unlikely (0.05% chance) that any field of or- 
der lO'xlO' includes the progenitor of a moderately rich 
cluster (volume density p — 10~^Mpc~'^), even in the 
long pencil beam from z —IS of a radio-SMG survey. 
This apparent paradox can only be resolved by reduc- 
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Fig. 2. — (top panel) Rcdshift histograms in GOODS-N (bins of 
5z = 0.02). from the surveys of Reddy et al. (2006) in UV-selected 
BX/BM galaxies, Team Keck Redshift Survey (Wirth et al. 2004), 
and Hawaii redshift survey (Cowie et al. 2004), the SMG surveys of 
Chapman et al. (2003, 2005), and the SFRG surveys of Chapman 
et al. (2004) and Casey et al. (2008). The z=1.99 structure can be 
seen as the highest peak of galaxies in the z ~ 2 region. The selec- 
tion functions for SMGs and BX/BMs are depicted normalized to 
the number of sources in each sample. The SFRGs are included in 
the SMG selection function (see text). There are a few overlapping 
BX/BM, TKRS and SMGs/SFRGs which are effectively double- 
counted in the plot. Note that the BX and SMG/OFRG galaxy 
samples in the z = 1.99 structure are completely non-overlapping, 
and the structure could arguably be more significant than other BX 
galaxy peaks in the combined BX/SMG galaxy sample, (lower 
panel) A zoomed in view of the structure with smaller redshift 
bins (Sz = 0.005). The spike galaxies listed in Table 1 are shown 
by the vertical bars. 

ing the frequency of such associations by a factor of ten 
(p = 10~^Mpc~'^ observed for moderately rich clusters 
compared with 10~^Mpc~^ implied by SMGs). Blain et 
al. concluded that either the dark matter halo masses 
are less extreme than inferred from a picture of sim- 
ple 1-1 biasing at a given mass, or some other biasing 
effect must be at work. Adelberger et al. (2005a) pre- 
sented an alternative explanation in the inherent bias of 
the cluster calculation methodology, however statistical 
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analyses of the SMG clustering through other methods is 
consistent with the Blain et al. (2004) calculation (e.g., 
Blake et al. 2006). The alternative explanation, that 
SMGs often reside in less overdense environments (and 
that ULIRG clustering could be dominated by the co- 
ordination of highly active periods across modest mass 
environments), is tested here, using UV observations of 
the richest overdensity found in Blain et al. (2004) and 
Chapman et al. (2005). 

In §2 we describe the sample properties, with §2.1 as- 
sessing the various selection functions, and a developing 
a statistically robust measure of the galaxy overdensity. 
In §2.2 we analyze the properties of the galaxies lying 
in the overdense structure, including their star forma- 
tion rates at a number of wavelengths, their luminosity 
function, and their angular positions on the sky. In §3 
we discuss the theoretical expectations of the measured 
redshift space contrasts, propose an explanation for the 
differential bias seen at different wavelengths (§3.1), and 
discuss the implications for future surveys at submm 
wavelengths (§3.2). We assume a A-CDM cosmology 
(Spergel et al. 2003) with flo = 0.3, flA ^ 0.7 and 
Hq = 70kms~^ Mpc""'^, so that 1 arcsec corresponds to 
8.4 kpc physical size at z = 2.0. 

2. CHARACTERIZING THE STRUCTURE AT UV- AND 
SUBMM- WAVELENGTHS 

Blain et al. (2004) highlighted an association of five 
SMGs in the GOODS-N field, aU lying within 1200 km/s 
of z = 1.995, representing the strongest association in 
their spectroscopic survey of 73 SMGs. The GOODS- 
N field is one of the better sampled fields from the ex- 
panded C05 SMG spectroscopic survey and the ongo- 
ing spectroscopic surveys of the SHADES submm fields 
(e.g., Coppin et al. 2008, Blain et al. in prep), how- 
ever associations of similar significance can be reasonably 
ruled out in these other survey fields (although the as- 
sociation distribution function N(iismg) is not yet well 
characterized). To further study this GOODS-N field 
SMG association and its properties at other wavelengths, 
we bring together an expanded sample of SMGs in this 
field, representing our primary sample of ultra-luminous 
star forming galaxies at high redshift. A population of 
submm- faint, radio-selected galaxies (SFRGs) with sim- 
ilar inferred star formation rates to SMGs as evidenced 
by radio and 24/im luminosities (Chapman et al. 2004, 
Casey et al. in prep.) represent a secondary sample from 
which to study the cluster and test our hypotheses. A 
map of source positions is shown in Fig. 1, depicting the 
spatial distributions of all galaxies with known spectro- 
scopic redshift z > 1 in the field and highlighting the 
candidate z = 1.99 structure. 

The GOODS-N sample of 30 SMGs comprises the spec- 
troscopic redshift sample of 22 radio-identified SMGs in 
Chapman et al. (2005), and six additional radio-SMGs 
with spectroscopic redshifts presented in Pope et al. 
(2006, 2008), and two further radio-SMGs from Daddi 
et al. (2008). We also include 19 submm-faint SFRGs 
from Chapman et al. (2004) and Casey et al. (in prep.), 
and assemble a complete catalog of UV/optical selected 
galaxies with redshifts, combining the surveys of Wirth et 
al. (2004), Cowie et al. (2004), Reddy et al. (2006), with 
our own spectroscopic data. In Table 1 we list all known 
galaxies in the z = 1.99 structure, including the several 
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Fig. 3. — Upper windows: Comparison of the Reddy et al. 
(2008) infrared (IR) luminosity function for z = 1.9 — 2.7 (trans- 
lated into SFR using Kennicutt 1998) to the SFRs for SMGs (red) 
derived from 24^m, and for UV-selected galaxies derived from ex- 
ponentially declining star formation models (blue). The IR lumi- 
nosity function of just the UV-selected galaxies in GOODS-N is 
also shown (not corrected for completeness) as squares, for direct 
comparison to the spike galaxies (which also are not corrected for 
completeness at the low luminosity end). As expected, the vast ma- 
jority of the overdensity compared to the average field appears at 
the high-SFR end in SMGs. Below is shown the raw distributions 
in SFR for UV-galaxies and SMGs/SFRGs. Lower windows: 
For the z~ 1.99 spike galaxies, distributions in radio (1.4 GHz) 
flux and 24/im flux density are presented, to show consistent mea- 
sures across UV-selected and SMG/SFRG populations, unaffected 
by uncertainties in SFRs calculated by different methods. 

unpublished UV-selected galaxies from our surveys. We 
assemble the relevant archival multi-wavelength data in 
Table 1 which is discussed explicitly in this paper, as well 
as new 1.4 GHz VLA radio measurements from the maps 
of Morrison et al. (2008). The redshift distributions are 
shown in Fig. 2. The well studied GOODS-N region is 
covered by a superb multi-wavelength survey to great 
sensitivity in most wavebands, and so it is unlikely that 
further missing galaxy populations could skew our results 
significantly. However there is no spectroscopic sample 
available which effectively targets quiescent galaxies at 
z - 2. 

2.1. Significance of the structure 

There is an obvious z ~ 1.99 galaxy concentration 
in the SMG redshift distribution lying in the GOODS- 



4 



N field which coincides with an apparently less signifi- 
cant concentration in the UV-selected galaxies. Figure 1 
shows the angular distribution of the galaxies, reveal- 
ing that the 9 SMGs and SFRGs in a redshift spike are 
not obviously clustered spatially across the field shown 
in Figure 1, spanning a region 7' (3.5 Mpc) on a side, al- 
though there are two relatively close pairs (of 9" and 17", 
76 kpc and 143 kpc respectively). The redshift distribu- 
tions around the z ~ 2 spike are compared in Figure 2. 

We first aim to characterize this structure with respect 
to galaxies selected at both UV and submm/radio wave- 
lengths. For the SMG sample, we define a selection func- 
tion by smoothing the redshift distribution of the ~100 
SMGs in the expanded survey from COS described above. 
We adopt a cluster finding algorithm based on an Erlan- 
gian distribution function (Eadie et al. 1971) 

p{Az\N\) = \{\Azf'^exp{-\Az)/iN-2y. (1) 

which characterizes the probability (in the absence of 
clustering) that a group of N galaxies would span a red- 
shift interval Az, considering that the product of our 
selection function and the density of SMGs is A. The 
procedure is very similar to that described in Steidel et 
al. (1998), and was used for our studies of halo substruc- 
tures in the Andromeda galaxy (e.g.. Chapman et al. 
2006). 

First looking just at the six SMGs in the structure, 
we find a probability of 0.005% for the concentration of 
SMGs at < z >= 1.992 with Az = 0.008 to be drawn 
from the selection function by chance, which maximizes 
the contrast with a peak of 5 of the 6 SMGs lying within 
the region. We define a galaxy overdensity as the num- 
ber of galaxies found in the structure (excluding the two 
galaxies which define the edges of the redshift interval) 
divided by the expected number in this redshift interval 
from the selection function (which is is only 0.4), yield- 
ing Sg = 4/0.4 = 10. This result is obviously dependent 
on the precise redshift adopted, as some scatter is intro- 
duced when individual galaxies are defined by different 
combinations of lines from molecular, neutral and ion- 
ized gas, which can differ by as much as 1000 km/s in 
a given SMG (Chapman et al. 2005, Greve et al. 2005, 
Swinbank et al. 2004). If the SMGs and SFRGs are con- 
sidered together, the contrast is maximized with 7 of the 
9 galaxies, and the probability to find the structure by 
chance rises marginally to 0.007%, with an implied over- 
density Sg = 8. 

Second, for the UV-selected galaxies, we adopt the se- 
lection function defined by the smoothed distribution of 
all galaxies in the spectroscopic survey of Steidel et al. 
(2004) and Reddy et al. (2006): see Fig. 1. The prob- 
ability of the z ~ 1.99 concentration being drawn at 
random from the selection function is 11.3%, and this 
contrast is maximized with a peak of 15 galaxies in the 
< z >— 1.987, Az = 0.015 region, with an overdensity 
Sg — nz/{n — 1) = 2.5. While C05 showed that many 
SMGs do have the colours of UV-selected galaxies, these 
spike SMGs and SFRGs do not overlap at all with the 
UV-selected galaxies in this z ~ 1.99 structure. 

We have therefore identified a highly significant over- 
density in ultra-luminous SMGs which appears not to 
correspond to a similarly overdense structure in less lumi- 
nous, UV-selected galaxies. The UV-galaxies have much 
lower SFRs and longer lifetimes (Reddy et al. 2006) than 
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Fig. 4. — Top windows: Comparison of histograms of stellar 
masses for SMGs in the spike, and for UV galaxies derived from 
continuous star formation models of Reddy et al. (2006), both in 
and out of the spike. The median values and bootstrap errors are 
shown with lines. Neither the UV galaxies nor the fSMGs show 
increased stellar mass buildup compared with the field, although 
the shear presence of the 9 SMG/SFRGs in the spike suggests a 
higher volume averaged stellar mass at z = 1.99 than in the field. 
Bottom window: For the z~ 1.99 spike galaxies, distributions 
in 4.5/im AB-mag (IRAC channel 2), representing a rest-frame 
~1.6^m or H-band light, showing the same trend as the stellar 
mass calculated above from the complete SEDs. 



SMGs, and it is a reasonable hypothesis that the UV- 
population would represent a more dependable tracer 
of large-scale structure. While there remains the pos- 
sibility that redshift incompleteness in quiescent galaxy 
populations (for example) could be responsible for the 
discrepancy, we will proceed with our explicit findings 
of a different bias between the two star forming galaxy 
populations. 

2.2. Analysis of the z = 1.99 spike galaxies 

The environmental context of high-redshift galaxy 
clusters is a topic which has lately received interest, since 
the preferential buildup of stellar mass in clusters should 
happen at earlier epochs than for field galaxies. Stei- 
del et al. (2005) found considerable excess stellar mass 
buildup and longer mean stellar ages in the densest struc- 
ture from their z ~ 2 surveys, as compared with the 
coeval field population. Clustering results from large 
spectroscopic surveys (e.g., VVDS -Meneux et al. 2008), 
rest- if mass-selected surveys using Spitzer-IRAC (e.g., 
de La Torre et al. 2007), and red galaxy selection (e.g,. 
BzK/ERO galaxies - Daddi et al. 2003) have come to 
similar conclusions regarding these environmental depen- 
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dencies. Our naive explanation of the differential bias 
between the subnim and UV-selected galaxies invokes a 
modest mass cluster in an active formation period, and 
we might therefore expect to find little in the way of 
stellar mass or age differences in the spike population. 
However, our interpretation could be clouded by missing 
populations of evolved galaxies which are not well sam- 
pled by GOODS-N spectroscopic surveys. We therefore 
analyze the spike galaxies to see if they show environ- 
mental effects characteristic of other high-^ clusters. 

First, we look at SFRs for all the galaxies in the spike 
(Fig. 3). Since many of the UV-selected galaxies are not 
detected in the 24/im or radio bands, while the SMGs 
and SFRGs typically have their SFRs underestimated 
by large factors in the UV, we calculate SFRs for the 
UV-galaxies from the dust-corrected UV luminosity, and 
from the 24/j,m luminosity for the SMGs and SFRGs. For 
the spike galaxies we also show the distributions of radio 
(1.4 GHz) and 24/im flux density to demonstrate that 
our derived SFRs are consistent with the raw data. The 
clustered environment allows for a useful comparison of 
SFR indicators without K-correction and model galaxy 
uncertainties. Our new radio measurements for the spike 
galaxies detect significant numbers of the lower lumi- 
nosity UV-selected galaxies in addition to all the SMGs 
and SFRGs, and radio-derived SFRs are generally con- 
sistent with those derived from Spitzer 24/im fluxes. We 
then consider the luminosity function (LF) of the whole 
z ~ 1.99 spike and compare with the LF derived for the 
whole UV population by Reddy et al. (2006) in Fig. 3. 
As expected for a relatively insignificant structure in the 
UV, the Reddy et al. LF is reasonably consistent with 
the z ~ 2 field galaxies, with a small excess expected 
from the overdensity factor of order two. When added to 
the spike sample, the SMGs clearly form a significant ex- 
cess compared with the average LF of all (dust-corrected) 
UV-selected galaxies. 

Next, we consider the stellar masses of galaxies in the 
spike (Fig. 4), compared to all non-spike UV-selected 
galaxies. The UV-galaxies in the structure do not exhibit 
the enhanced buildup of stellar mass which is apparent 
in other high-z UV/optical-selected clusters (Steidel et 
al. 2005, Kurk et al. 2008). The UV spike galaxies have 
<M*>= 1.1±3.6 X lO^" Mq, while the "field" UV galax- 
ies have <M*>= 1.5 ± 3.4 x 10^° Mq (derived from the 
tables in Reddy et al. 2006). Using the stellar masses 
estimated for the SMGs by Hainline (2008), Hainline et 
al. (in prep.), we find that the spike SMGs clearly have 
larger stellar masses (<M*>= 6.0 ±3.6 x 10^" Mq) than 
the UV-selected galaxies, but the SMGs' stellar masses 
are not significantly different from the median for 73 
SMGs from Hainhne (2008) (M* = 6.9 x 10^° Mq). We 
note that Borys et al. (2005) derived even larger (5x) 
stellar masses for these spike SMGs. Hainline (2008) as- 
cribe the difference in their masses to those from Borys 
et al. as due to two issues. Firstly the use of different 
stellar evolution models: the Maraston TP-AGB models 
(Maraston et al. 2006), and differing light- to-mass and 
SF histories, and secondly the removal of AGN contribu- 
tions to the observed 8/im flux prior to fitting. Alexan- 
der et al. (2008) have also explicitly attempted to correct 
for the AGN contribution to the rest-frame K-band light 
in these SMGs and derive lower masses (factor of ^ 2). 
We stress that the relative masses of the BX/BMs and 



SMGs will be robust against some of the systematic un- 
certainties in the models, and that our results are not 
sensitive to changing the masses by factors of a few, and 
hence we believe our conclusions are robust. We also 
show in Figure 4 the 4.5/im fluxes of the spike galax- 
ies (rest-frame ~ 1.6/xm or i?-band) revealing the SMGs 
are clearly much brighter than the UV-galaxies with- 
out uncertainties in stellar mass model fitting. There 
is therefore evidence for significant stellar buildup in our 
z ~ 1.99 overdense environment, although this is driven 
exclusively by the large overdensity of SMGs. 

Comparing our "spike" SMGs to "typical" SMGs at 
similar epochs (e.g.. Chapman et al. 2005, Borys et al. 
2005, Pope et al. 2006, Hainline et al. 2008) reveals no 
obvious differences in SFR or stellar mass (the median 
values of both are similar for spike and field SMGs - 
<SFR>=810 Me/yr, <M* >=7xl0i" Mq). An addi- 
tional route to studying the activity in the structure ver- 
sus the field would be through AGN populations. Cata- 
loged AGN in GOODS-N (e.g., Cowie et al. 2005, Reddy 
et al. 2007) do not reveal any significant enhancement at 
z ^ 2, nor do the Xray properties of these spike SMGs 
(Alexander et al. 2005). 

3. DISCUSSION 

3.1. Evolutionary context of the structure 

To place our SMG cluster in a cosmological context, we 
must assume a galaxy bias for SMGs. We adopt h ^ 2.1, 
a value consistent with simulations in van Kampcn et 
al. (2005), and similar to the estimated galaxy bias in 
the UV galaxy sample (Adelberger et al. 2005b). Fol- 
lowing the method introduced by Steidel et al. (1998), 
we relate the true mass overdensity 5 to the observed 
redshift-space overdensity of SMGs Sg through the ex- 
pression 

l + b6=\C\il+5^g) (2) 

where |C| = Vapp/Vtme takes into account the effects of 
redshift-space distortions from peculiar velocities. For 
our SMGs, SFRGs and UV galaxies, Vapp ^ 9900 Mpc^ 
is the co-moving volume in which the measurement is 
being made, bounded by the lO'O by lO'O region on the 
plane of the sky and the co-moving distance, neglecting 
peculiar velocities, between z = 1.978 and z = 2.008. 
The volume correction factor can be approximated by 
C = 1 + f - f{l + (5)1/3 where / = fi„j(z)°'5 = 0.96 at 
z ~ 2.0. Thus S ~ 2.3 with = 10.0, and C = 0.53. The 
linear overdensity, ^o, can then be approximated using 
equation 18 (spherical collapse) from Mo & White (1996) 

<5o ^l(l + (5)-'/' + ^2(l + <5)"°-''''^ + A3(l + (5)-l/2+A4 

with Al = -1.35, A2 = 0.78785, ^3 = -1.12431, A4 = 
1.68647. Thus for our z - 1.99 SMG structiu^e {S = 
2.3), we find 60 ~ 0.85 as the real space linear matter 
overdensity. 

We now address what the SMG overdensity will have 
become by redshift z ^ 0. Evolved forward in our 
adopted cosmology, its linear over-density of Sq ^ 0.85 
at z = 2 corresponds to a linear over-density of (5p ~ 2.13 
at z = 0. This exceeds the collapse threshold 5c = 1.69. 

However, if we repeat the analysis with the UV- 
selected galaxy Sfjy ^ 2.5 case, we find a matter over- 
density at z = 2 of only 5 = 0.83, a linear overdensity 
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do = 0.5 or Sp ^ 1.1 at z = 0, and the structure asso- 
ciated with this apparently modest overdensity of UV- 
galaxies would not have collapsed by 2; = 0. 

The mass of a virialized object, as we would infer for 
the SMG overdensity, can be calculated by assuming that 
the SMG overdensity is associated with an Eulerian mat- 
ter overdensity of 6m — 2.3 at z = 2.0 (e.g., Steidel et 
al. 1998), so its mass is (1 + Sm)pVtrue ^ 10^^ Mq where 
p is the mean co-moving matter density and Vtme = 
Vapp/C ~ 18700 Mpc^ is approximately the volume we 
calculated within the observed overdensity. There is, 
however, no straightforward way to estimate the mass 
of the lower contrast UV-galaxy spike, as it has not ob- 
viously virialized by z = 0. Wc simplistically assume the 
inferred mass would be >3x lower. 

3.2. Explaining the differential bias 

Our task is then to explain why such a significant struc- 
ture in the SMG population exhibits such a modest over- 
density in the UV-selected galaxies. One possibility is 
that the true overdensity and cluster mass concentration 
is in fact as large as implied by the SMGs, and we have 
missed a large population of evolved quiescent galaxies, 
selected for instance through the restframe 2/xm emis- 
sion. At a redshift z ~ 2, such a population is difficult 
to confirm spectroscopically, and might explain for in- 
stance why we found no difference in the stellar masses 
of field versus spike UV-selected star forming galaxies. 
However, wc sec no logical reason why the mass concen- 
tration should then be represented by SMGs (and not 
UV-galaxies) , which are an even more luminous popula- 
tion of star forming galaxies. 

Another line of inquiry is to look for mechanisms to 
produce the apparent differential bias directly. Dynam- 
ical masses have been estimated for galaxies selected in 
the restframe UV using nebular emission line widths 
(Pettini et al. 2000; Erb et al. 2006), and for SMGs 
and SFRGs using both nebular lines and rotational CO 
molecular lines (e.g., Fraycr et al. 1998; Neri et al. 2003; 
Swinbank et al. 2004, 2006, 2008; Chapman et al. 2008). 
Attempts to account for the fact that estimates of the 
dynamical masses of high-z galaxies are significantly less 
than the clustcring-infered halo masses (Adelberger et 
al. 2005b; Blain et al. 2004) have fallen under the names 
of temporal bias (Scannapieco & Thacker), merger bias 
(e.g., Furlanetto & Kamionkowski 2006) or assembly bias 
(Dalai et al. 2008) . These approaches attempt to account 
for the rapid mergers and other properties of galaxies in 
halos boosting the apparent clustering signal, and they 
might naturally extend to the scenario we have uncov- 
ered in this z ~ 2 structure, which represents an even 
more radical version of this problem. As pointed out by 
Furlanetto & Kamionkowski (2006), the extended Press- 
Schechter and Mo & White (1996) biasing scheme yields 
no merger bias as it explicitly ignores the variation of 
merger rates with the large-scale density field. However, 
for populations like SMGs, even an analytic approxima- 
tion to the bias implicit in aggressively merging galax- 
ies is suggestive that the clustering length coiild be en- 
hanced. The challenge is to apply this formalism not just 
to a pair- wise enhancement of the population, but to the 
synchronization of 9 SMG-class galaxies spread over an 
angular region 3.5 Mpc on a side. 

We consider from Furlanetto & Kamionkowski (2006), 



equation 24, which describes how pairs of peaks in the 
density field are clustered, as a function of the clustering 
of the peaks themselves (i.e., identifying the peaks with 
galaxies, and the "pairs" as recently merged objects). In 
the limit that the peaks are well above the characteristic 
mass scale of objects (likely a reasonable assumption for 
SMGs), the correlation function of pairs x(^) is 

X(r-) = [1+C,(r)r-1 (3) 

where £,g(r) is the correlation function of the galaxy 
population of interest {bg is the associated matter cor- 
relation function). This gives an enhancement to the 
overall clustering of the underlying galaxy population. 

Ignoring redshift space corrections, we then make the 
following argument: there is one true matter overdensity, 
Sm- Suppose that UV-galaxies and SMGs are drawn from 
halos with the same mass (and hence bias buv)- Then, if 
SMGs correspond to pairs of objects, they have an extra 
"pair" bias bp = ^/(x/^g)- So the observed overdensities 
5" will be: 

1 + Sljy = 3.5 = 1 + buvSm (4) 
l + S^SMG = ^^ = ^ + bpbuv5m (5) 

which requires 6j, ~ 4 (slightly different if redshift space 
correction is included) . 
While the model is approximate and the observations 

are not necessarily complete, the result is not unrea- 
sonable. A remaining question is what is the correct 
r which applies in eq. 2&3 above to evaluate bp. On very 
large scales, where « 1, we can Taylor expand eq. 3 
about r = 0, x(r) ~ 4^g(r), to infer a pair bias bp ~ 2. 
Our physical interpretation of this is two pairs (or four 
SMGs), so we obtain four factors of £,g. However, on 
smaller scales, the pair bias goes up because one is deal- 
ing with peaks: on scales where the correlation function 
of SMGs is large, ^ 1 (i-c, where clusters are form- 
ing), x(^) = 15, so bp ^ 4. This is similar to our required 
factor of boosting, and likely sufficient given uncertain- 
ties in the observations. 

We emphasize that simulations that have tried to quan- 
tify this "assembly bias" give smaller effects, ~ 10% 
(Wechsler et al. 2002; Zao et al. 2003a,b; Dalai et al. 
2008). However, these studies haven't looked at this ex- 
act question of recent mergers, which may be a stronger 
effect than the concentration, age, spin, etc., which 
are modeled in these simulations, and therefore our toy 
model is not inconsistent with these studies. 

3.3. Context and Implications for SCUBA2 surveys 

While the most overdense forming environments at 
z > 1.5 have been observed to have large numbers 
of ultra-luminous galaxies associated (e.g., the SSA22 
z = 3.09 proto-cluster has 8 known SMGs - Geach et al. 
2006), our results suggest that if one searches for galaxy 
clusters using the most vigorously bursting (and typically 
highest luminosity) systems, one is prone to finding the 
brief periods when systems of lower matter overdensity 
are rapidly evolving. The exceptionally strong concen- 
tration of SMGs presented here corresponds to a modest, 
not statistically significant (in terms of 5) concentration 
in UV-selected galaxies, whereas the expectation would 
be that this submm "spike" should represent one of the 
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largest mass structures known. While we have outlined 
a toy model to explain the scenario in § 3.2, it remains 
to be seen how such scenarios can be realized by detailed 
theoretical models and simulations of cluster formation. 

Baugh et al. (2005) have presented a model of the SMG 
population in a semi-analytical galaxy evolution frame- 
work. Using the the Millennium Simulation (Springcl et 
al. 2005), Almeida et al. (2008 in prep.) predict that the 
correlation function of model SMGs from Baugh et al. 
(2005) is close to a power law {rj = (r/ro)'') over more 
than three decades in separation, with 7 = —1.94 ± 0.05 
and a correlation length ro = 8.8 ± 0.3 Mpc, slightly 
smaller but consistent with that measured by Blain et al. 
(2004) from the redshifts of SMGs, ro = 9.8 ± 3.0 Mpc. 
Swinbank et al. (2008) performed a detailed compari- 
son of actual SMG observations with model SMGs in 
the Baugh et al. (2005) model. In the Millennium Sim- 
ulation, this correlation length corresponds to a pre- 
dicted halo mass of Mhaio = ^-^tH x 10^^ Mq (Gao 
et al. 2005). Measuring the dark halo mass for the 
model SMGs directly, Swinbank et al. (2008) derive 
Mhaio = 3.6^f:5 X 10^^ Mq. Taken together this suggests 
there is no merger bias in the correlation lengths: the 
correlation length of the SMGs is identical, within the 
errors, with the whole population of halos of the same 
masses. Our discovery in this paper of a large submm 
burst in a modest mass environment suggests that a more 
thorough investigation is needed of the clustering prop- 
erties of SMGs using combined semi-analytic and numer- 
ical simulations to trace the star formation and merger 
histories of high-redshift galaxies. 

We note that similar criticisms could apply to the clus- 
tering of AGN (e.g., Gilli et al. 2005), and overdensi- 
ties selected through AGN (e.g., Rottgering et al. 1996, 
2003; Brand ct al. 2003). In particular, the overdensi- 
ties of other galaxies around luminous AGNs arc typ- 
ically poorly characterized relative to the field. Our 
well calibrated results on SMGs suggest that these AGN- 
identified galaxy clusters could be far less massive than 
suggested. 

While our characterization of a single cluster of SMGs 
cannot be extrapolated to cluster formation in general, 
it is noteworthy that large scale surveys in the submm 
with SCUBA2, and the Herschel and Planck space tele- 
scopes, and even existing wide-field surveys at 24/xm from 
Spitzer (e.g., Farrah et al. 2006) may reveal rather mod- 
est mass structures associated with some of the highest 
contrast source overdensities. The large beam size of 
Planck and the ffersc/ie/-Spectral and Photometric Imag- 
ing REceiver (SPIRE) 500/um survey may be particu- 
larly sensitive to this sort of effect (Negrello et al. 2005, 
2007), where the low-resolution surveys will measure the 
summed contributions of groups of sources within the 
beam. Further, the number of these low-mass but active 
clusters are likely significant enough to make it impossi- 
ble to associate a large ro with a large dark matter halo 
mass. Detailed foUowup of the highest concentrations 
of ultra-luminous sources, although arduous, is clearly 
warranted. In the future, serendipitous searches for CO 
emitters at the same redshift using ALMA, and evidence 
for boosted flux from confused sources using Herschel 
may also provide insight into the nature of this struc- 



ture. 

4. CONCLUSIONS 

We have described an apparent cluster of SMGs and 
RGs at < z 1.99 in the GOODS-N field, the strongest 
known association of SMGs (Blain et al. 2004; Chapman 
et al. 2005). The implied matter over-density from the 
redshift space contrast Sg ~ 10 is is ~ 10^^ M©. Search- 
ing the same volume in UV-selectcd galaxies, we find a 
much reduced overdensity ( 6g ~ 2.5). The low SFRs of 
UV-galaxies and much longer star formation timescales 
on average compared with the SMGs and SFRGs sug- 
gest they should be more dependable for defining the 
mass scale of this structure. Wc conclude that we are 
observing a modest mass overdensity (which likely will 
not even virialize by z = 0), which is undergoing a par- 
ticularly active star formation phase. 

We have characterized the bolometric luminosity func- 
tion of the z ^ 1.99 spike galaxies compared to that 
estimated for a large sample of UV-selected star form- 
ing galaxies (Reddy et al. 2008), finding that it is almost 
entirely in the SMG population where an excess is seen. 
The clustered environment allows for a useful comparison 
of SFR and stellar mass indicators without K-correction 
and model galaxy uncertainties. Our new radio measure- 
ments for the spike galaxies detect significant numbers 
of the lower luminosity UV-selected galaxies in addition 
to all the SMGs and SFRGs, and radio-derived SFRs 
are generally consistent with those derived from Spitzer 
24/xm fluxes. 

While the most overdcnsc forming environments at 
z > 1.5 have been observed to have large numbers of 
ultra-luminous galaxies associated, our results here sug- 
gest that if one searches for galaxy clusters using the 
highest luminosity systems, one is prone to flnding the 
brief periods when systems of lower matter overdensity 
are rapidly evolving. SMGs don't obviously trace the 
most massive structures in the Universe, as simulations 
and constraints to date have suggested they should. The 
success of our toy model with merger bias explaining the 
basic elements of our observations suggests that simula- 
tions may not currently be accurately reflecting the pro- 
cesses occurring in the formative phases of galaxy clus- 
ters. Further work to embed these aspects of merger 
bias in simulations will therefore be of great interest in 
understanding galaxy cluster assembly. Investigation of 
similar systems in upcoming surveys at long wavelengths 
will be of great interest to explore the full range of en- 
vironments associated with strong overdensities of ultra- 
luminous galaxies. 
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TABLE 1 

Properties of all UV-selected galaxies, SMGs and SFRGs in GOODS-N z = 1.99 spike" 



ID 


RA 


dec 


rcdshift'' 


type of rcdshift 


Sl.4GHz 


S24jim 


MAB,5.8f4m 


MAB,4.5Mm 


bMG-yj 


12 35 00.13 


+62 10 47.2 


1.994 


UV, Ha 


128.5+8.1 (26.5) 


1237+8.0 


nn oo 1 n n^ 

20.32+0.27 


nn n^ I n no 

20.97+0.08 


bMG-132 


12 36 18.32 


+62 15 50.5 


1.999 


CO(4-3) 


172.0+8.4 (34.3) 


344.0+8.0 


20.31+0.18 


20.65+0.07 


SMG-140e 


12 36 21.25 


+62 17 08.3 


1.994 


Spitzer 


169.4+8.8 (31.9) 


370.0+11.4 


20.06+0.10 


20.31+0.06 


SMG-140W 


12 36 20.99 


+62 17 09.5 


1.989/1.988 


UV,Ha 


63.4+10.6 (9.0) 


d 


20.67+0.09 


20.96+0.07 


SMG-169 


12 36 32.53 


+62 07 59.7 


1.993 


UV 


80.4+8.6 (15.1) 


820.8+15.2 


19.64+0.15 


20.29+0.07 


bMG-172'^ 


12 36 35.57 


+62 14 24.0 


2.001 


UV 


77.0+7.8 (16.3) 


1410.0+13.9 


18.29+0.09 


18.90+0.06 


bMCi-255 


12 37 11.99 


+62 13 25.6 


1.992 


UV, CO 


50.2+8.1 (10.2) 


nor' n I ^ r\ 

225.0+7.0 


m '^n 1 n m 

21.79+0.21 


n "1 A K 1 n nn 

21.45+0.09 


SFRG-130 


12 36 17.54 


+62 15 40.7 


1.993 


UV 


234.3+8.6 (45.1) 


20.2+5.2 


20.77+0.17 


20.58+0.08 


SFRG-179 


12 36 40.73 


+62 10 11.0 


1.968/1.977 


UV 


72.5+8.3 (13.7) 


83.0+7.8 


21.29+0.19 


21.33+0.07 


SFRG-254 


12 37 11.32 


+62 13 30.9 


1.996/1.993 


Ha/CO(4-3) 


127.4+8.7 (23.0) 


537.3+9.1 


19.72+0.11 


20.01+0.06 


BX-1201 


12 36 14.13 


+62 15 41.8 


.../2.000 


UV 


10.6+7.0 (2.8) 


29.0+5.5 


23.04+0.27 


22.99+0.07 


BX-1192 


12 36 16.83 


+62 15 14.3 


.../1.996 


UV 


21.3+10.3 (3.1) 


111.6+10.7 


21.47+0.15 


21.56+0.07 


BX-1071 


12 36 20.10 


+62 11 12.6 


.../1.996 


UV 


< 31.8 






23.06+0.21 


BX-1228 


12 36 20.19 


+62 15 40.6 


1.999/1.995 


UV 


< 14.1 


14.7+4.1 


23.05+0.22 


23.10+0.11 


BX-1267 


12 36 22.67 


+62 16 21.6 


1.996/1.996 


UV 


< 14.5 


48.7+7.2 






BX-1339 


12 36 25.09 


+62 17 56.8 


1.993/1.984 


UV 


< 14.9 




22.63+0.31 


23.03+0.12 


BX-1307 


12 36 48.33 


+62 14 16.7 


.../2.002 


UV 


18.2+6.4 (5.2) 


145.5+12.2 


21.34+0.09 


21.61+0.07 


BX-1329 


12 36 54.62 


+62 14 07.7 


.../1. 987 


UV 


< 18.2 








BX-1827 


12 36 56.63 


+62 15 19.0 


1.988/... 


UV 


9.6+4.1 (2.3) 






23.46+0.13 


BX1129 


12 36 56.94 


+62 08 48.7 


.../1.973 


UV 


21.7+11.2 (3.0) 


104.4+10.3 


21.46+0.09 


21.60+0.07 


BX-1431 


12 36 58.48 


+62 16 45.5 


2.006/1.996 


UV 


< 13.4 






23.25+0.34 


BX1378 


12 37 02.02 


+62 14 43.4 


.../1.971 


UV 


10.4+4.2 (2.5) 


29.4+5.5 


22.43+0.10 


22.40 0.07 


BM-1122 


12 37 02.62 


+62 11 56.7 


1.994/1.986 


UV 


< 17.6 






23.30+0.16 


BX-1434 


12 37 16.80 


+62 14 38.8 


.../1.994 


UV 


4.6+6.1 (1.5) 






23.18+0.27 


BM-69 


12 37 1 7 68 


+62 14 359 


1 991 /I 991 


UV 


< 17.8 


QO 7-1- K O 






BX-184 


12 37 19.28 


+62 13 00.6 


1.998/... 


UV 


< 12.2 






22.95+0.25 


BM-70 


12 37 20.05 


+62 14 57.1 


1.997/1.994 


UV 


< 30.9 






23.69+0.19 


BX-1605 


12 37 21.51 


+62 18 30.6 


1.977/1.970 


UV 


< 23.1 






23.36+0.17 


BM-1155 


12 37 23.98 


+62 12 12.1 


2.024/2.015 


UV 


< 33.6 






23.28+0.13 


BX-1642 


12 37 32.40 


+62 17 50.8 


2.010/2.004 


UV 


< 21.0 






22.95+0.09 


BX-1708 


12 37 32.65 


+62 19 10.6 


.../1.987 


UV 


16.3+9.1 (3.0) 


208.3+14.6 


21.03+0.07 


21.23+0.07 


BX-1694 


12 37 33.82 


+62 18 46.3 


2.009/2.005 


UV 


10.1+8.4 (2.0) 


13.8+4.3 


22.28+0.11 


22.58+0.07 


SCBX-6106 


12 37 10.82 


+62 08 19.0 


.../1.984 


UV 


< 13.5 




22.41+0.55 


23.15+0.22 


SCBX-6946 


12 36 57.70 


+62 10 36.8 


.../1.992 


UV 


< 14.7 




22.53+0.31 


22.84+0.25 



" Galaxies listed have redshifts which lie within ±2a of the z = 1.99 peak, and as shown by the vertical bars in Fig.2b. 

Redshifts of UV-selected galaxies are quoted as emission line/absorption line, where ... indicates that one or the other was not 
possible from the spectrum - see Reddy et al. (2006) for details of the BX/BM galaxies, whereas we have adopted this same procedure 
for the SCBX galaxies (from our own redshift surveys) , selected from a color-cut transformed to AB-mag from that presented in 
Steidel et al. (2004). Molecular gas (CO) redshifts are from Chapman et al. (2008), and Small et al. (in prep). The Spitzer-IRS 
redshift comes from Pope et al. (2008). 

^ Radio fluxes from Morrison et al. (2008), derived from the full resolution (1.7") map. Integrated flux measurements and errors are 
shown, with the signal-to-noise in brackets. For UV-galaxies, source fluxes are listed if a radio peak lies within 0.5" of the optical 
source, otherwise a 3cr upper limit. 

24^m flux density of SMG-H140w is confused with SMG-H140e, although the radio emission is clearly resolved by MERLIN 
(Chapman et al. 2004), and even by the VLA (Biggs & Ivison 2006). 

HDF172 could conceivably be interpreted as an SFRG rather than an SMG. It is a SCUBA photometry detection in Chapman et al. 
(2001), a SCUBA-mapping non-detection in Borys et al. (2004), and a weak (~ 3a) MAMBO-mapping source (at 1.2mm) from 
Greve et al. (2008). 



